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Abstract 

In 60 pb~^ of data taken on the ^-'(3770) resonance with the CLEO-c detector, we find 8 — > 
fi^v event candidates that are mostly signal, containing only 1 estimated background. Using this 
statistically compelling sample, we measure preliminary values of B{D^ M^^) = (3-5 ± 1.4 it 
0.6) X 10-^, and determine /^+ = (201 ± 41 ± 17) MeV. 



*Submitted to the 32"'^ International Conference on High Energy Physics, Aug 2004, Beijing 
^On leave of absence from University of Chicago. 
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I. INTRODUCTION 



Measuring purely leptonic decays of heavy mesons allows the determination of meson 
decay constants, which connect measured quantities, such as the BB mixing ratio, to CKM 
matrix elements. Currently, it is not possible to determine experimentally from leptonic 
B decays, so theoretical calculations of /b must be used. The most promising of these 
calculations involves lattice QCD [1] [2] [3], though there are other methods [4] [5] [6] [7] [8]. 

Measurements of pseudoscalar decay constants such as //)+ provide checks on these cal- 
culations and help discriminate among different models. 




FIG. 1: The decay diagram for /i+z/. 
The decay diagram for fi^iy is shown in Fig. 1. The decay rate is given by [9] 

r(z)+-/+-) = ^/WMz..(i-^)Vc/ , (1) 

where Md+ is the mass, mi is the mass of the final state lepton, Vcd is a CKM matrix 
element equal to 0.2205 [10], and Gp is the Fermi coupling constant. Various theoretical 
predictions of fD+ range from 190 MeV to 350 MeV. Because of helicity suppression, the 
electron mode D'^ — > e'^u has a very small rate. The relative widths are 3.2 : 1 : 2.4 x 10~^ 
for the T~^iy, jjC^v and e^v final states, respectively. Unfortunately the mode with the largest 
branching fraction, t^v, has at least two neutrinos in the final state and is difficult to detect. 

II. THE CLEO-C DETECTOR 

The CLEO-c detector is equipped to measure the momenta and direction of charged 
particles, identify charged hadrons, detect photons and measure with good precision their 
directions and energies. Muons above 1.1 GeV can also be identified. The detector is almost 
cylindrically symmetric with everything but the muon detector inside a superconducting 
magnet coil run at a current that produces an almost uniform 1.0 T field. The detector 
consists of a six-layer wire drift chamber at small radius that is low mass, suitable for these 
relatively low energies. It is followed by a 47-layer drift chamber; both chambers use a gas 
mixture of 60% Helium and 40% Propane. These two devices measure charged track three- 
momenta with excellent accuracy. The drift chamber also measures energy loss, dE/dx, that 
is used to identify charged tracks below about 0.7 GeV/c. After the drift chamber there is a 
Ring Imaging Cherenkov Detector (RICH) [11] , that identifies charged particles over most 
of their momentum range. The RICH is surrounded by a Thallium doped Csl crystal array 
consisting of about 8000 tapered crystals 30 cm long and about 5x5 cm^ at the rear. 
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III. DATA SAMPLE AND SIGNAL SELECTION 



In this study we use 60 pb~^ of CLEO-c data recorded at the resonance (3.770 GeV). 
These events consist of a mixture of D'^D~ , D°D and continuum events. There also may 
be small amounts of r+r" pairs and two-photon events. 

We examine all the recorded events and retain those containing at least one charged D 
candidate in the modes listed in Table I. The selection criteria are described in detail in 
what follows. 

All acceptable track candidates must have a helical trajectory that approaches the event 
origin within a distance of 0.005 m in the azimuthal projection and 0.05 m in the polar 
projection, where the azimuthal projection is in the bend view of the solenoidal magnet. 
Each track must posses at least 50% of the hits expected to be on a track and its polar 
direction must have cosine with respect to the beam direction of magnitude < 0.93. 

We use both charged particle ionization loss in the drift chamber (dE/dx) and RICH 
information to identify kaons and pions used to fully reconstruct D mesons. The RICH is 
used for momenta larger than 0.55 GeV. Information on the angle of detected Cherenkov 
photons is translated into a Liklihood of a given photon being due to a particular particle. 
Contributions from all photons associated with a particular track are then summed to form 
an overall Liklihood denoted as Ci for each particle hypothesis. To differentiate between 
pion and kaon candidates, we use the difference: —2 log(£7r) + 2 log(£ft:). Usually this cut is 
set at zero except for muon candidates where the difference between —2 log(i2^) + 2 log(£x) 
is set at 10, to ensure a high, well understood efficiency. To utihze the dE/dx information 
we calculate (7^, as the difference between the expected ionization loss for a pion and the 
measured loss divided by the measurement error. Similarly, Gk is defined in the same manner 
using the expected ionization for a kaon . 

We use both the RICH and dE/dx information for D~ meson tag candidate tracks in the 
following manner: (a) If neither the RICH nor dE/dx information is available, then the track 
is accepted as both a pion and a kaon candidate, (b) If dE/dx is available and RICH is not 
then we insist that pion candidates have PIDdE = o'n ~ '^k < kaon candidates have 
PIDfiE > 0. (c) If RICH information is available and dE/dx is not available, then we require 
that PIDrich = -2 \og{C^) + 2 \og{CK) < for pions and PIDrjch > for kaons. (d) If 
both dE/dx and RICH information are available, we require that {PIDdE + PI D rich) < 
for pions and (PIDdE + PIDrich) > for kaons. 

We reconstruct 7r°'s by first selecting photon candidates from energy deposits in the 
crystals not matched to charged tracks that have deposition patterns consistent with that 
expected for photons. Pairs of photon candidates are kinematically fit to the nominal 7r° 
mass. We require the pull, the difference between the raw and fit mass normalized by its 
uncertainty, to be less than three for accepted 7r° candidates. 

K° are formed from a pair of charged pions which are kinematically fitted to come from 
a single vertex. We also require that the invariant mass of the two pions be within 4.5 times 
the width of the Kg mass peak which is 0.004 GeV r.m.s. 

IV. RECONSTRUCTION OF CHARGED D TAGGING MODES 

Tagging modes are fully reconstructed by first evaluating the difference in the energy of 
the decay products with the beam energy. We then normally require the absolute value of 
this difference to be within 0.02 GeV of zero, approximately twice the r.m.s. width, and 
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then look at the reconstructed D beam constrained mass defined as 

mo = ^^£Lr^"(E^> (2) 

where i runs over all the final state particles. We also use the charge-conjugate tags and 
search for — > iJi~v^] in the rest of this paper we will not mention the charge-conjugate 
modes explicitly, but they are always used. 

The tud distribution for all tagging modes considered in this data sample are shown in 
Fig. 2 and listed in Table I along with the number of signal and background events. The event 
numbers are determined from fits of the tjid distributions to Gaussian signal functions plus 
a background shape parameterized as 3'''^ order polynomial for the K^'k~'k~'k°, KsTI^ti^-k^ 
and Ks7r~7i° and from fits to double Gaussian signal functions plus S^'^ order polynomial for 
K~^n~n~ and K^tt" tags where we see a small tail on the higher mass side. 



Mode 


Signal 


Background 


K~^Tr~TT~ 


15188 ± 233 


583 ± 336 




4082 ± 81 


1826 ± 343 


KsTT- 


2110 ±81 


251 ± 19 


K s'IT~ TT'^ 


3975 ± 81 


1880 ± 342 




3297 ± 87 


4226 ± 516 


Sum 


28652 ± 285 


8765 ± 784 



TABLE I: Tagging modes and numbers of signal and background events determined from the fits 
shown in Fig.2. 

Selecting those candidates within 3 r.m.s. widths of the D~ mass reduces the signal 
number by 77 events giving a total of 28575±286 events used for further analysis. 

V. D+ n+Uf, SELECTION CRITERIA 

To select l^^^ii events we first reconstruct event candidates and then search 

for events with a single additional charged track presumed to be a /u^. Then we infer the 
existence of the neutrino by requiring a measured value near zero (the neutrino mass) of the 
missing mass squared (MM^) defined as 

MM^ = {E,,am - E^+f - (-to- - t^^+f , (3) 

where Jf^- is the three-momentum of the fully reconstructed D~ . 

We need to restrict the sample to candidate /i"*"//^ events resulting from the other D. 
Thus we wish to exclude events with more than one additional charged track, which we 
take to be the muon candidate, or with extra neutral energy. It is possible, in fact even 
hkely, that the decay products of the tagging D~ interact in the detector material, mostly 
the EM calorimeter and spray tracks and neutral energy back into the rest of the detector. 
To evaluate the size of these contributions we use a very pure sample of events obtained 
by finding fully reconstructed D°D events. The numbers of these events in various decay 
modes are fisted in Table II, we use a total of 782 events. 
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FIG. 2: Beam Constrained Mass distributions for different fully reconstructed D~ decay candi- 
dates; the curves show the sum of Gaussian signal functions and S^'^ order polynomial background 
functions. A single signal Gaussian is used for all modes except for (a) and (c) where two Gaussians 
are used, (a) D~ A'+tt^tt", (b) i^+vr^vr^vr", (c) KsTt^, (d) KsTT~tt^tt~^ 

and (e) D~ i^svr^vr'^. (preliminary) 
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Mode 1 Mode 2 # of events 







89 




~Tr~ K~Tr~^ 


392 


K+7^-7^^ 


~Tr~ K~ TT'^ TT~ -K~^ 


301 



TABLE II: Fully reconstructed D°D° events 

The number of interactions of particles with material and their consequences depend 
on the number of particles, the kind of particles and their momenta. Thus, the sum over 
these neutral D decay modes isn't quite the same as the sum over the tagging D~ decays, 
however the average over these modes is quite similar to the D~ tagging modes for this level 

of statistics. 

Extra tracks do appear in these D"D events. None of these tracks, however, approach 
the main event vertex. Requiring that good tracks are within 0.05 m along the beam and 
0.005 m perpendicular to the beam does not include any additional tracks from interactions 
in the material. D~ tags with additional n~^7r~ candidates are rejected. 

In the D°D events, energy in the calorimeter not matched to any of the charged tracks is 
shown in Fig. 3. Figure 3(a) shows the energy of the largest shower and 3(b) shows the total. 
We accept only as extra showers those that do not match a charged track within a connected 
region. A connected region is a group of adjacent crystals with energy depositions which 
are nearest neighbors. This suppresses a lot of hadronic shower fragments which would 
otherwise show up as unmatched showers. Hadronic interactions and very energetic tt^'s 
tend to produce one connected region with many clusters. For further analysis we require 
that the largest unmatched shower not to be larger than 250 MeV. This requirement is 
(93.5±0.9)% efficient for signal events. 

The muon candidate is required to be within the barrel region of the detector | cos 6 \ < 
0.81, where 6 is the angle of the muon with respect to the beam electrons; this requirement 
insures that the MM^ resolution is good as tracks at larger angles with respect to the beam 
are measured with poorer precision. In addition, this requirement helps reject background 
from the decay tt+tt"; this mode also gives a MM^ near zero. Requiring the muon 

candidate in the barrel region (the tt"*" in this case) avoids having the photons from this 
decay being lost in the transition region of the calorimeter between the barrel and the 
endcap, because the tt" direction is almost directly opposite the tt"*". Furthermore, the muon 
candidate is required not to be consistent with the kaon hypothesis using RICH information. 
Finally, we also require that the muon candidate deposits less than 300 MeV of energy in 
the calorimeter, characteristic of a minimum ionizing particle. This requirement is very 
efficient for real muons, and rejects about 40% of the pions as determined using a sample 
of reconstructed D° K^n^ decays. Fig. 4 shows the muon deposited energy in the EM 
calorimeter both from data on e^e^ — > fji^ n~ and from GEANT simulation of the same 
process. The Monte Carlo and data are in excellent agreement. We therefore use a GEANT 
simulation of — > /i+i/ with lower energy muons to determine that the efficiency of the 
calorimeter energy cut is (98.7 ± 0.19)%. 

Table. HI shows the properties of each muon candidate from the 8 events in the signal 
region. When evaluating MM^ using equation 3 there are two important considerations 
that are not obvious. First of all, we exphcitly need to take into account the crossing angle 
between the and e" beams. This angle is about 4 mrad, varying slightly run to run; we 
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FIG. 3: Largest and Total Extra shower energies in the D°D° sample. 
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FIG. 4: Muon deposited energy in CC of Muon pairs from (a) DATA and (b) MC. 
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Tag 


MM2 


CC energy 


-21og(£i^) 


-21og(£^) 




(GeV2) 


of /i+(GeV) 








0.032 


0.186 


-4.3 


-166.0 + 


KsTT 


-0.019 


0.201 


0.00 


-140 


KlTTT 


-0.051 


0.190 


31.9 


-252.9 + 


Ktttt 


-0.004 


0.221 


0.00 


-115.2 + 


KsTTTt" 


0.032 


0.164 


-0.32 


-130.6 - 


KsTTTTTT 


0.001 


0.245 


-11.7 


-138.9 + 


Kinnr" 


0.002 


0.204 


-8.6 


-88.6 




0.014 


0.208 


-8.3 


-113.0 + 



TABLE III: Muon Candidate Properties. 

use this information and Lorentz transform all laboratory quantities to the center-of-mass. 
Secondly, we require that the reconstructed D~ have exactly the known D~ mass. This 
changes and improves somewhat our knowledge of the direction. 

The MM^ from Monte Carlo simulation is shown for our different tagging samples in 
Fig. 5. The signal is fit to a sum of two Gaussians with the wider Gaussian having about 
30% of the area independent of tagging mode. The resolution [a) is defined as 

a = /lai + (1 - /i)a2, (4) 

where ai and o"2 are the individual widths of the two Gaussians and fi is the fractional area 
of the first Gaussian. The resolution is approximately 0.025 GeV^ consistent among all the 
tagging decay modes. 

We check our simulations using the D' — > Kgi^" decay. Here we choose events with 
the same requirements as used to search for //'•"z/ but require one additional found Kg. The 
MM^ distribution for this final state is shown in Fig. 6. The resolution is measured to be 
0.024±0.002 GeV^ from a single Gaussian fit, consistent with but slightly larger than the 
Monte Carlo estimate of 0.021±0.001 GeV^. To account for the difference in resolution 
between data and simulations we scale the resolution by 14 % to 0.028 GeV^ when looking 
for the — > /J^^'^ signal. 

The MM^ distributions for our tagged events requiring no extra charged tracks besides 
the muon candidate and showers above 250 MeV as described above is shown in Fig. 7. We 
see a small signal near zero containing 8 events within an interval that is twice as wide as our 
resolution, -0.056 GeV^ to -1-0.056 GeV^. This signal is most likely due to the fJ~^^iJ. 
mode we are seeking. The large peak centered near 0.25 GeV^ is from the decay K°7r~^ 
that is far from our signal region. 

A typical event is shown in Fig. 8. 

VI. BACKGROUND EVALUATION 
A. Introduction 

There are several background sources we need to evaluate. These include background from 
other modes, background from misidentified D°D events and continuum background. 
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FIG. 5: Monte Carlo simulation of l^^T^fi events for different tags. The plots have been fitted 

to two Gaussians centered at zero where the second Gaussian constitutes around 30% of area. 

The requirement of the muon depositing <300 MeV in the calorimeter, while about 99% 
efficient on muons, rejects only about 40% of pions as determined from the D°D event 
sample where the pion from the K^ti"^ mode was examined. In Fig. 9 we show the deposited 
energy in the calorimeter for both kaons and pions obtained from the Kti tag sample. 

B. Backgrounds 

There are a few decay modes that could simulate the signal. These are listed in 
Table IV along with the background estimate we obtained by Monte Carlo generation and 



10 



30 
20 
10 



0.00 0.10 0.20 0.30 0.40 0.50 

MM2(GeV2) 



160 
120 
80 
40 



O.OO 0.10 0.20 0.30 0.40 0.50 
l\/ll\/l2(GeV2) 



FIG. 6: MM^ distribution for the decay Ksir' from data and MC 

reconstruction of each specific mode. The branching ratios are from the Particle Data Group 
except for the 7r"'"7r° mode where a separate CLEO analysis gives a somewhat lower value 
[12]. This mode is the most difficult to reject because the MM^ peaks very close to zero, at 
0.018 GeV^, well within our resolution of 0.028 GeV^. While we have insisted that the muon 
candidate be well within our acceptance, it is possible for the photons from the n° decay to 
inadvertently be matched to the tracks from the tagging D~ or be missed. The maximum 
photon energy of the n" at the generator level is shown in Fig. 10. We note that at least 
one photon from the tt+tt" mode exceeds our 250 MeV calorimeter energy requirement and 
should in most cases cause such a decay to be vetoed. 

Even though the i^^yr"*" mode gives a large peak in the MM^spectrum near 0.025 GeV^, 
our simulation shows that only a very small amount can enter our signal region, only 0.06 
events. We have simulated backgrounds from t^v. Out of 10,000 simulated events 

with D~ tags, we found background only when r+ n^u. Because of the small D~^-t~^ mass 
difference, the r"*" is almost at rest in the laboratory frame and thus the vr"*" has relatively 
large momentum causing the MM^ distribution to populate only the low MM^ region, even 
in this case with two missing neutrinos. The MM^ distribution is shown in Fig. 11. 

The semileptonic mode 7r°/i+z/^ is similar to tt^tt" except that the 7r° often carries off 
enough momentum to result in large MM^. We found no candidate background events in a 
Monte Carlo sample consisting of 50,000 tags plus a —>■ Ti^ji^v decay. 
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FIG. 7: MM^ using D~ tags and one additional opposite sign charged track and no extra energetic 
showers (see text). The insert shows the signal region for — > /x+z/ enlarged, ±2(7 range is shown 
between the two arrows, (preliminary) 

Mode B (%) # of events 
tt+tt" 0.13±0.02 0.31±0.04 
K°7r+ 2.77±0.18 0.06±0.05 
T+iy 3.2x n+iy 0.36±0.08 
TT^n+u 0.31 ± 0.15 negligible 



TABLE IV: Backgrounds from specific decay modes 

C. D°D° and Continuum Backgrounds 

These backgrounds arc evaluated by analyzing Monte Carlo samples corresponding to 15 
times the total amount of data in our possession. To normalize our Monte Carlo events 
to our data sample we used cr^o^" ~ 2. 9 nb and (Tcontinuum = 6.4 nb. In each sample we 
found one background event within two standard deviations of zero. These correspond to 
0.16±0.16 D°D events and 0.17±0.17 continuum events forming background. 

D. Background Summary 

Our total background is 1.07±0.25 events. Because of the uncertainties in the Monte 
Carlo simulation we assign a 100% error to our background estimate: 1.07±1.07 events. 

VII. BRANCHING RATIO AND DECAY CONSTANT 



Subtracting the 1.07 event background from our 8 events in the signal region, we deter- 
mine a branching fraction using a detection efficiency for the single muon of 69.9%. This 
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Run: 202742 Event: 98595 

Ks7i~jf'"7i^ Tag 




FIG. 8: Typical D ^ u^^ event. The muon is shown as well as the two charged pions forming 
the Kg. We notice the presence of a curler track (a pion) with momentum around 50 MeV. 

efficiency includes the tracking, the particle identification, the probability of the crystal en- 
ergy being less than 300 MeV and the probability of not having another unmatched shower 
in the event with energy greater than 250 MeV. We assign a relative 5.3% error on this 
efficiency, the components of which are shown in Table V. We use a 3% systematic error 
on track finding found using the double tagged events and we estimated the error on the 
particle identification cut to be 1% from studies of D*^ decays in higher beam energy data. 
The error on the minimum ionization cut is evaluated from the efficiency of the cut by 
generating signal Monte Carlo. To evaluate the error on the extra shower cut we evaluated 
the efficiency of that cut in the D°D sample discussed above. The total systematic error is 
obtained by summing all entries in quadrature. 

To compute the branching ratio we use 6.93 signal events divided by 69.9% and the 
28575±286 tags. No other efficiencies enter. We find 

B{D+ /i+z/^) = (3.5 ± 1.4 ± 0.6) x 10"^ , (5) 

The error on the background contributes 15.4% to the systematic error on the branching 
ratio. 
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FIG. 9: Deposited energy in EM calorimeter for (a) Kaons, (b) Pions from D' 





Systematic error % 



MC statistics 0.8 

Track finding 3 
PID cut 1 
Minimum ionization cut 1 

Extra showers cut 4 
Number of tags 1 



TABLE V: Systematic errors on —> jJi^fn efficiency. 

The decay constant /£>+ is then obtained from Equation 1 as 

/^+ = (201 ± 41 ± 17) MeV . 
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FIG. 10: Maximum photon energy of the tt" in the D tt tt" decay. 



VIII. CONCLUSIONS 



There have been several experimental studies of D meson decay constants. The Mark III 
group published an upper limit of B{D'^ f^^^fi) < 7.2 x 10~^, which leads to an upper 
hmit on the decay constant /d+ < 290 MeV at 90% confidence level based on 9.3 pb~^ of 
data taken on the i/j" [13]. BES claimed the observation of one event at a center-of-mass 
energy of 4.03 GeV with a branching ratio of (0.08lo;o5)*^ [15]. Recently, using 17.7 pb~^ 
of tp" data they presented a 3 event signal and claimed a background of 0.25 events where 
neither n'^ir", or t'^v were mentioned as a possible background mode, nor was continuum 
background considered [16]. Here they find a branching ratio of (0.12to;og3lo;oo9)%, and a 
corresponding value of f]j+ = {'iG^lllllllll) MeV. Our value is considerably smaller, though 
compatible with their large error. 

Our analysis shows the first statistically compelling signal for — > //+z/. The prelimi- 
nary branching fraction is 

B{D+ n+v) = (3.5 ± 1.4 ± 0.6) x 10"^ , (7) 

and the preliminary decay constant is 

= (201±41± 17) MeV . (8) 

Our result for at the current level of precision, is consistent with predictions of all 
of the models listed in Table VI. 
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FIG. 11: Missing Mass squared distribution for t^u and r"*" — > ir^v. 



Model 



/^+ (MeV) 
210 ± lO+jl 



fD+/fD+ 

1.13 ± 0.02+^ "^ 



Lattice (unquenched) (UKQCD) [3] ^j.^^ _i_ j.w_jg _i_ u.u^_q q2 
Lattice (partially quenched) (MILC) [2] 215 ± el^^t^l^ 1.14 ± O.Ollg:g^ ± 0.03 ± 0.01 

QCD Spectral Sum Rules [5] 203 ± 20 1.15 ± 0.04 

QCD Sum Rules [6] 195 ± 20 

Relativistic Quark Model [7] 243 ±25 1.10 

Potential Model [4] 238 1.01 

Isospin Mass Splittings [8] 262 ± 29 



TABLE VI: Theoretical predictions of fD+ and f^+/fjj+ 



The models generally predict /^+ to be 10-15% larger than /£>+. CLEO previously 
measured /^+ as (280 ± 19 ±28 ±34) MeV [17], and we are consistent with these predictions 
as well. We look forward to more data to improve the precision. 
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